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ABSTRACT 

This paper examines the ultraviolet and X-ray photons generated by hot gas in the Galactic thick 
disk or halo in the Draco region of the northern hemisphere. Our analysis uses the intensities from 
four ions, C IV, O VI, O VII, and O VIII, sampling temperatures of ~ 10 5 to ~ 3 x 10 6 K. We measured 
the O VI, O VII, and O VIII intensities from FUSE and XMM-Newton data and subtracted off the 
local contributions in order to deduce the thick disk/halo contributions. These were supplemented 
with published C IV intensity and O VI column density measurements. Our estimate of the thermal 
pressure in the O Vl-rich thick disk/halo gas, pth/k = 6500t 2 60u K cm ~ 3 , suggests that the thick 
disk/halo is more highly pressurized than would be expected from theoretical analyses. The ratios 
of C IV to O VI to O VII to O VIII intensities were compared with those predicted by theoretical 
models. Gas which was heated to 3 x 10 6 K then allowed to cool radiatively cannot produce enough 
C IV or O Vl-generated photons per O VII or O VHI-generated photon. Producing enough C IV and 
O VI emission requires heating additional gas to 10 5 < T < 10 6 K. However, shock heating, which 
provides heating across this temperature range, overproduces O VI relative to the others. Obtaining 
the observed mix may require a combination of several processes, including some amount of shock 
heating, heat conduction, and mixing, as well as radiative cooling of very hot gas. 
Subject headings: Galaxy: general — Galaxy: halo — ISM: general — ISM: individual (Draco) - 
ultraviolet: ISM — ultraviolet — X-rays 



1. INTRODUCTION 

Observations of C IV, O VI, O VII, and O VIII 
show that hot gas extends a few kpc ab ove the plane 
(|Savage fe Wak kcr 200§ lYao et al.l 120091 into a region 
historically called the halo and more recently called 
the thick disk. This is an active and dynamically im- 
portant region of the Galaxy. On-going star forma- 
tion has be en seen in the thic k disks of other spi- 
ral galaxies (jReuff fc Howkll2010h . implying that these 
thick disks continue to experience the winds and su- 
pernova (SN) explosions that result from young stars. 
In addition, thick disks churn as ou tflows from the 
thin disk, long postulated theoretically ([Shapiro fe Field! 
11976ft and more recently se en in numerical simula- 
tions having normal SN rates dJoung fe Mac Low! [20061 : 
Ide Avillez fe Breitschwerdtl 120071) . make their way into 
the thick disk. Estimates based on the observed C IV, 
O VI, 1/4 keV, and 3/4 keV intensities at high southern 
latitudes indicate a large radiative loss rate; the region 
above the thin disk is so luminous that it radiates 2/ 3 
of the Galaxy's supernova power ([Shelton et al.1 12007T ). 
Since most of the Galaxy's SN energy is injected nearer 
to the midplane, a large luminosity produced above the 
thin disk suggests the sort of churning seen in simula- 
tions. Not only is the thick disk affected by in situ heat- 
ing and upwellings from the thin disk beneath it, but it is 
also affected by inflows from above, as in the case of high 
veloc ity clouds colliding with the Galaxy ([Tripp et al.l 
2003), or if we consider the more extended halo, in the 
form of accretion from the intergalactic medium (e.g., 
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iRasmussen et al.1l2009h . 

Most of the volume within a few kpc of the midplane is 
filled with very hot gas having tem peratures in excess o f 
10 6 K (see pathlength esti mates in IShelton et al.l (|2007f ) 
or scale height estimates in lYao fc Wangil2007l ). This gas 
is traced by O VII ions and O VIII ions having collisional 
ionization equilibrium (CIE) temperatures, Tcie, of ~ 
1 x 10 6 K and ~ 3 x 10 6 K, respectively. Somewhat cooler, 
but still hot, gas is traced by C IV ions {T C ie ~ lx 10 5 K) 
and O VI ions {Tqie ~ 3 x 10 5 K) . It fills a lesser fra ction 
of the space (see pathlengths in Sh elton et al.ll2007l ) but 
accounts for more radiative energy loss than does the 
hotter gas. Because gas in this temperature regime cools 
rapidly, it must be replenished from a nearby reservoir of 
hotter gas. For this reason, a common conception of the 
hot interstellar medium (ISM) is one in which 1 to 3 x 
10 5 K gas resides in transition z ones between hotter and 
cooler gas. Some authors (e.g.. ISavage fc Wakkerl[2009l l 
have begun to call the 1 to 3 x 10 5 K material "transition 
temperature" gas. Here, for simplicity, however, we use 
the term "hot gas" for the entire ~ lx 10 5 to ~ 3 x 10 6 K 
regime. 

The tracers of hot gas are observed by UV and X- 
ray instruments. Ultraviolet instruments are used to ob- 
serve the strong resonance line transitions (2s 2 Siy 2 - 2p 
2 P 3/2 and 2s 2 Si /2 - 2p 2 P 1/2 ) of C IV and O VI. Ow- 
ing to the excellent spectral resolution of recent UV in- 
struments, interstellar C IV and O VI have been studied 
well via absorption line spectroscopy (ISavage et aL1l2003l : 
iBowen et all 120081 ISavage fc Wakkerl 120091) Thev have 
also been seen via emission spectroscopy, but along fewer 
sight lines (e.g IShelton et al.1 l200lt [Dixon fc Sankritl 
2008; Pa rk et all [20091. In contrast, X-ray instruments 
are used to observe Ka transitions of O VII and O VIII. 
Due to the limited spectral resolution and/or through- 
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put of current X-ray instruments, it is far easier to de- 
tect photons produced by interstellar O VII and O VIII 
ions than to identify O VII and O VIII absorption 
profiles in the spectra of distant X-ray point sources. 
While som e absorption column densities have been mea- 
sured (e.g.. lYao fe W ang 20071: iBregman fc Lloyd-Daviesl 
120071) . most of the O VII and O VIII measurem ents are of 
emitted photons (e.g.. iHenlev fc Sheltonll2010l ). 

Column density and emission intensity measurements 
plumb different aspects of the interstellar gas. The col- 
umn density depends only on the path length and the 
density of pertinent ions in the gas. The emission in- 
tensity, however, depends on two additional factors, the 
electron density and the emissivity, which is a function 
of temperature. Thus, high ions that are much cooler 
than their CIE temperature (for example, ions result- 
ing from photoionization and ions that have cooled in 
a non-equilibrium fashion) can be seen via absorption 
spectroscopy, but are very inefficient radiators. This is 
especially relevant for C IV, because a significant frac- 
tion of the C IV ions seen in absorption measurements 
are photoionized ([Savage & Wakkcr 2009). 

Here, we focus on the emission intensities from C IV, 
O VI, O VII, and O VIII from a single region of the sky, 
(£ = 90°, b ~ 40°), a region near, but not toward, the 
Draco molecular clo ud complex (MBM 41 through 44, 
lMagnani~e t al. 1985). This paper is organized as fol- 
lows. In Section [21 we assemble measurements of C IV, 
O VI, O VII, and O VIII. The C IV intensity measurements 
presented in Subsection 12.11 come from published obser- 
vations made by FIMS/SPEAR. We complement these 
with new measurements of the O VI intensity, made from 
FUSE data, (Subsection |2~2|) . In Subsection [2~3l we de- 
duce the extraplanar O VI column density in the Draco 
region from existing FUSE measurements. In Subsec- 
tion [23] we present our O VII and O VIII measurements, 
made from archival XMM-Newton observations of a di- 
rection 1° from our FUSE and FIMS/SPEAR direction. 
For each ion, we take care to remove the local contri- 
butions (from the Local Bubble and/or heliosphere) and 
compensate for interstellar absorption. The results are 
the second known foursome of intrinsic C IV, O VI, O VII, 
and O VIII intensities for the region of the thick disk/halo 
that is above the thin disk and toward a single part of the 
sky, and the first such foursome for the northern Galactic 
hemisphere. In Section [31 we use the O VI emission in- 
tensity and absorption column density to determine the 
density and pressure of the O Vl-rich gas in the thick 
disk/halo above the thin disk. We then use the ratios of 
our C IV, O VI, O VII, and O VIII intensities to test phe- 
nomenological models. These include a radiative cooling 
model, such as would be expected for accretion of in- 
tergalactic gas onto the galaxy, a point-injection model, 
such as would be expected for supernova explosions, and 
a turbulent mixing model. We find that radiative cool- 
ing of very hot gas is not able to produce enough C IV 
and O VI photons per O VII or O VIII photon. Bubbles 
formed from point injections of 5 x 10 50 ergs have the 
opposite problem. Models of turbulent mixing between 
very hot gas and cooler gas do not provide O VII or O VIII 
predictions, but can explain high ratios of C IV to O VI 
emission. We also present the ratio of C IV, O VI, O VII, 
and O VIII intensities, so that they may be compared 
with other model predictions as they become available. 



Table 1 

Observed and Extraplanar Intensities of C IV, O VI, O VII, and 



O VIII 


Ion 


i,b 


Observed 


Intrinsic Extraplanar 






Intensity 


Intensity 




Degrees 


(ph s _1 cm" 2 sr — 1 ) 


(ph s~ 1 cm" 2 sr" 1 ) 


C IV 
O VI 
O VII 
O VIII 


90.0°, 39.6° 
90.0°, 39.6° 
90.0°, 38.4° 
90.0°, 38.4° 


5500 ± 1830 
3500 ± 940 
fi s+ - 9 

i o+ - 4 


6540 ± 2180 
4770S 

1.0 ± 1.2 



Note. — The extraplanar O VI column density in this region is 
N OVI = 1.9 ± 0.5 x 10 14 cm" 2 . 

The C IV and O VI intensities are for the doublets. 

In Section [31 we summarize our results, concluding that 
energy must have been injected into the halo in a way 
that heated some of the gas to ~ 3 x 10 6 K while heating 
other gas to 1 to 3 x 10 5 K. 

2. DATA 
2.1. The Halo C IV Intensity 

As part of a survey prog ram, the Draco region was ob- 
served by FIMS/SPEAR ([Edelstein et al.l I2006T ). which 
is sensitive to C IV resonance line doublet photons 
(A = 1548 1 551 A ), between February and May of 2004. 
iPark et all (|2009l ) presented measurements of the dou- 
blet's intensity as a function of location within the Draco 
neighborhood in the form of a smoothed, pixelated map. 
The direction of interest, I = 90.0°, b = 39.6° is located 
near the intersection of four 0.1°-wide pixels in their map. 
For these four pixels, the average intensity is ~ 5500 in 
units of ph s -1 cm -2 sr _1 , hereafter line unit s (LU), and 
the ra nge of intensities is 700 LU. Because IPark et all 
( 2009) reported a signal to noise ratio of > 3.0 for inten- 
sities over 5000 LU, we take la to be 1/3 of the measured 
intensity and adopt / = 5500 ± 1830 LU as the observed 
intensity in the doublet. This and later measurements 
are tabulated in Table [T] This value is ve r y nea r to the 
average intensity reported by IPark et all (|2009l ) for all 
off-cloud regions within their map: 5559 ± 1008 LU. 

No observations of foreground C IV emission have been 
reported in the literature and little emission is expected 
from the Local Bubble or heliosphere. Therefore, we 
assume that the C IV-rich gas resides above the Local 
Bubble, and also above the Galaxy's H I and/or dust- 
rich layer. The photons emitted in this region would 
be subject to extinction by intervening gas. The degree 
of obscuration determined from the H I column density 
from the Leid en-Argentine-Bonn ( LAB) survey of neu- 
tral hydrogen (|Kalberla et al.|[2005l ) is roughly similar to 
that determined from the DIRBE-corrected IRAS data 
(jSchlegel et al.l[l998h . The LAB survey found a column 
density of Nhi = 1-52 x 10 20 cm -2 , which implies a 
loss of 20% of the intensity, based on the color excess 
to H I column dens i ty rel ation empirically determined 
bv iDiplas fc Savage! (I1994T) and the extinction relation 
calculated bv iFitzpatrickl (|1999h . The DIRBE-corrected 
IRAS 100 fim intensity is 7i o = 0.71 MJy sr~\ which 
implies a loss of 1 1%, based on the Jig o to color excess 
conversion used bv lSchlegel et alj (119981 ) and the extinc- 
tion relation calculated bv lFitzpatrickT (fl999). Thus, if 
all of the obscuring material lies between the source of 
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the photons and the Earth, then the halo's intrinsic in- 
tensity would be 6200±2070 to 6880±2290 LU, assuming 
the 11% and 22% loss rate, respectively. Here, we take 
the halo's intrinsic C IV doublet intensity to be the av- 
erage: 6540 ± 2180 LU. Note that throughout the paper, 
the unrounded measurement values are used in the cal- 
culations, but the rounded results are reported. 

2.2. The Halo O VI Intensity 
2.2.1. FUSE Observations and Data Reduction 

We obtained two FUSE observations for this project. 
The first, observation ID El 2401, was directed toward 
the Draco cloud, t = 89.7°, b = 38.5°. Beginning on 
2005-11-30, six exposures were taken for a total of 5,144 
seconds. The FUSE spacecraft failed before the remain- 
der of the observation could be completed. The sec- 
ond observation, observation ID E12402, was directed 
toward an unobscured region near the Draco cloud, t = 
90.0°, b = 39.6°. Several exposures began on 2006-11-15 
and the remainder began on 2007-3-14. The total expo- 
sure time was 79,912 seconds, of which 31,316 seconds 
were taken at night. The original intent of this program 
was to perform a cloud shadowing study. However, the 
on-cloud observation was too short to yield clear detec- 
tions (see subsection 12.2. 2|) . Here, we focus on the off- 
cloud observation. 

We examined the data taken through the largest aper- 
ture of the LiFIA channel, as this yields spectra with 
the greatest signal to noise ratio. Like other observa- 
tions of diffuse emission, the data were taken in time-tag 
mode. Photon lists for each exposure were processed 
with version 3.2 of the CalFUSE pipeline. We used the 
extended source extraction window option, disabled the 
background subtraction feature, and set the pulse height 
limits to 2 and 25 for this processing. 

Earlier analyses found that the wavelength scale must 
be re-calibrated by comparisons with e mission lines of 
known wavelengths (jShelton et al.ll200ll ). Because each 
exposure may have a different offset, the exposures were 
processed by CalFUSE and re-calibrated before being 
added. The wavelength offsets were found by aligning the 
observed Lyman ft airglow feature in the LiFIA spectra 
with the theoretical Lyman (3 wavelength at zero velocity 
in the geocentric reference frame. After being corrected 
in this manner, the wavelength scale for each spectrum 
was converted to the heliospheric reference frame and the 
spectra added together. 

2.2.2. O VI Intensity Measurements 

Figure [1] shows the LiFIA spectrum in the 1027 to 
1040 A region. The O VI 1032 A resonance line is appar- 
ent in both the night-only and the day+night spectra. 
The O VI 1038 A resonance line is naturally dimmer and 
also easily obscured by the neig hboring C II* 1037 A fea- 
ture. The bright peaks near 1027.5, 1028, and 1039 A 
in the day+night spectrum are due to atmospheric hy- 
drogen and oxygen. The 1031 A feature that appears in 
the day+night spectrum, but not in the night-only spec- 
trum, is thought to be the second order diffraction line 
of atmospheric He I (|Shelton et al.ll2007t ). 

Here, we focus on the O VI 1032 A line in the off-cloud 
data. We measure its strength, both by the method of 
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Figure 1. LiF 1A spectra of the off-cloud observation. The spec- 
tra have been converted to the Local Standard of Rest reference 
frame and binned into 0.104 A pixels. The rest wavelengths of the 
O VI, C II*, 2nd order He I, and Earth airglow features (©) are 
marked on the plots in the LSR reference frame. Top: Data taken 
during the night portion of the satellite's orbit. Bottom: Data 
taken during both the night and day portions of the orbit. During 
the orbital day, the detectors are subjected to H I and O I airglow 
photons and also record second-order diffraction counts from He I, 
which appear around 1031 A. The rest wavelengths of the O VI 
and C II* emission lines are marked with vertical bars on the lower 
plot. 
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Figure 2. Method #2 fits to the spectra (left: night-only data; 
right: day+night data). The plotted spectra have been converted 
to the LSR reference frame and binned into 0.104 A pixels. The 
Method #2 fitting functions are overplotted. They consist of a 
linear function that models the continuum, a Gaussian convolved 
with a tophat that models the 1032 A feature, and in the case of 
the day+night data, a second Gaussian convolved with a tophat 
that models the airglow 1031 A feature. The ragged curves near 
the bottom of each plot are the error bars on the data. 

searching for significant exces ses in counts that was em- 
ployed in lShelton et al.l ()200lD and later papers (method 
#1) and by a modified version of the method of fitting 
with a Gaussian function that has been convolved with 
a 106 km sec" 1 to phat function that was employed in 
iDixon et all (|200H ) and later papers (method #2). The 
modification is to add a second tophat-convolved Gaus- 
sian to model the second order He I airglow feature that 
appears at 1031 A in the daytime data. See Figure [2] for 
the fitting results. Both methods cleanly separate the 
O VI 1032 A feature from the daytime 1031 A feature. 
Table [2] presents the intensity measurements. We aver- 
age the measurements, yielding an observed intensity in 
the O VI 1032 A line of 2340 ± 630 LU. 

To place the O VI 1032 A intensity along our off-cloud 
sight line into context, we compare it with O VI 1032 
A intensities measured for othe r high-latitude lin es of 
sight. Using arch i val FU SE data. IDixon et al.l (|2006| ) and 
IDixon fc Sankritl (|2008l ) have searched for diffuse O VI 
photons along roughly 300 lines of sight. Here we con- 
sider the 54 sight lines in their sample (plus E12402) 
with Galactic latitudes of \b\ > 20° and orbital-night 
exposure times exceeding 18 ksec. This sample is com- 
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Table 2 

Observed Intensities and 1 a Statistical Uncertainties for the 1032 A O VI feature 



Night Only 


Night Only 


Day+Night 


Day+Night 


Method #1 


Method #2 


Method #1 


Method #2 


(ph s _1 cm -2 sr — 1 ) 


(ph s -1 cm -2 sr -1 ) 


(ph s — 1 cm -2 sr -1 ) 


(ph s — 1 cm -2 sr - - 1 ) 


off-cloud: 1970 ± 600 


2240 ± 660 


2140 ± 470 


2990 ± 790 


on-cloud: 


± 690 


-240 ± 1600 


± 1030 



Note. — For the on-cloud data, method #1 was unable to fit the continuum for the night-only 
data. Method #2 calculated the 3<r upper limit. For comparison with the other values, the 3cr upper 
limits found by method #2 were converted into a "measurement ±1<t" form by dividing by 3. 



plete in the sense that any emission brighter than 2000 
LU would have been detected in an 18 ksec night-time 
exposure. Within this sample, 63% have statistically- 
significant (i.e., > 2a) intensities of O VI 1032 A pho- 
tons, the greatest of which is 5500 LU. If null detections 
are treated as having zero flux and a 1 a error, then the 
median and mean intensities in the 1032 A line are 2100 
and 1800 LU, respectively, while the standard deviation 
about the mean is 1620 LU. Our off-cloud sight line, with 
an observed 1032 A intensity of 2340±630 LU, lies within 
one sigma of the mean, so it is not at all unusual. 

Our measurements of the weaker O VI 1038 A feature 
were not as sound as our measurements of the stronger 
1032 A feature. The 1038 A feature sits to the right of 
the 1037 C II* line and to the left of the O I airglow fea- 
ture, making it impossible to determine the background 
continuum accurately. For this reason, we present nei- 
ther day+night nor night-only 1038 A measurements. 

We do, however, present measurements from our search 
for a second O VI 1032 A feature. This search was moti- 
vated by the observation of a weak intermediate velocity 
{Vlsr = — 112 km sec -1 ) feature in the LAB survey of 
galactic H I. If infalling gas at this LSR velocity was 
accompanied by O VI, then the stronger line in its dou- 
blet would appear at a wavelength of 1031.49 A in our 
heliocentric restframe spectra. We searched for such a 
feature in a 0.36 A wide window (corresponding to the 
width of the LWRS aperture) in the night-only data. We 
did not search for the feature in the day+night data, be- 
cause it would have overlapped with the second order 
He I feature. Using method #1, we found a statistically 
insignificant excess intensity of 210 ± 480 LU. 

We examined the on-cloud data, as well. These spec- 
tra were of very poor quality, because of the short obser- 
vation time (5144 seconds of day+night data, including 
2056 seconds of night-only data.) Neither method #1 
nor method #2 could find emission features around 1032 
or 1038 A. We were able to use method #1 to determine 
and subtract the continuum and then measure the resid- 
ual intensity within the wavelength range where the 1032 
A feature was expected to be found, based on the loca- 
tion of the 1032 A feature in the off-cloud spectra taken 
from the day+night data. This method yielded a mea- 
sured intensity of -243 ± 1604 LU, where the negative 
sign indicates that this region of the spectrum has fewer 
counts than the fitted continuum. Method #2 yielded 3 
a upper limits of 3090 and 2060 LU from the day+night 
and night-only data, respectively. In both cases, the un- 
certainties are too large to constrain our conceptions of 
the interstellar medium. 



Table 3 

Velocity with Respect to the Local Standard of Rest 



Night Only 


Night Only 


Day+Night 


Day+Night 


Method #1 


Method #2 


Method #1 


Method #2 


(km sec -1 ) 


(km sec -1 ) 


(km sec -1 ) 


(km sec" 1 ) 


O VI 1032 A -12 + 7 


2 + 17 


-9 + 11 


-3 + 14 



In order to determine the velocities using method #1, 
the heliocentric spectra were refit using a Gaussian for 
the feature in question and a low-order polynomial for 
the continuum. Method #2 also yielded velocity esti- 
mates. In both cases, the reference frames of the ana- 
lyzed spectra were heliospheric. The resulting velocities 
have been converted to the LSR reference frame by the 
addition of 16.73 km sec -1 . (See Table [3] for the mea- 
surement results.) The wavelength scale is accurate to 
- 10 km sec" 1 (|Shelton et al.ll200l . 

2.2.3. Intrinsic Halo O VI Intensity 

For optically thin plasmas, the 1038 A line is half as 
strong as the 1032 A line. Assuming this ratio, the dou- 
blet intensity becomes 3500 ±940 LU (Table [J). The ob- 
served intensity is due to photons emitted in the halo but 
attenuated by intervening material plus photons emitted 
locally, by the Local Bubble if it exists. When the obser- 
vations were planned, it was intended to use the on-cloud 
measurement to estimate the foreground intensity. How- 
ever, because the exposure time was too short for a good 
measurement, we take the local contribution from the 
only other O VI cloud shadow observation, that toward 
£ = 278.6°, b = -45.3° where a cloud lies 230±30 pc from 
Earth. That cloud observation yielded upper limits on 
the O VI 1032 A and 1038 A line intensities. The tightest 
upper limit is that on the 1032 A line in the day+night 
data. Multiplying it by 1.5 in order to account for the 
1038 A contribution yields a limit on the doublet in- 
tensity of 30±l*° LU (|Shelton|[200l . The unrounded 
values were made available by the author. Subtract- 
ing them from our unrounded doublet intensity yields 
3480lt« LU. 

If the non-local photons originated above even part of 
the Galaxy's H I layer, then the intrinsic intensity would 
be greater than 3480^0 

LU. Here we make the assump- 
tion that the O Vl-rich gas resides above the full extent 
of the Galaxy's H I and/or dust-rich layer. As with the 
C IV intensity, we estimate the degree of obscuration of 
the halo O VI emission using the LAB H I column density 
(jKalberla et al.l [2005T) and the DIRBE-corrected IRAS 
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100 fim intensity (jSchlegel et al.lH998l ). Using the same 
input values and conversion relations that we used in Sec- 
tion 2.1, we find that the LAB H I column density and 
the 100 /zm intensity imply O VI losses of 33% and 19%, 
respectively. Thus, if all of the obscuring material lies 
between the source of the photons and the Local Bub- 
ble, then the halo's intrinsic intensity would be 4320^2,™ 
to 5220ti5i° LU. Here, we take the halo's intrinsic O VI 
doublet intensity to be the average: 4770^ggQ LU, which 
is also listed in Tabled] 

The optical thickness of the on-cloud line of sight is 
much larger (Iioo = 3.48 MJy sr _1 , implying a loss of 
70% of the incoming photons). If the gas above the 
cloud produces the same doublet intensity as that de- 
rived above for the off-cloud line of sight, then only 
1440^420 LU would survive passage through the cloud 
and reach the detector. Of this, 960^280 LU would be 
in the O VI 1032 A line. This intensity is within our 
observational upper limits. 

2.3. An Estimate of the Halo's O VI Column Density 

We estimate the line of sight O VI column d ensity in 
nearly the same manner as IShelton et all (|2007| ). namely 
by averaging the column densities on the four nearest 
directions for which they are measured and calculating 
the uncertainty from the typical scatter in measured val- 
ues. The differenc e between our method and that of 
IShelton et al.l (|2007l ) is that we also include the measure- 
ment error; see below. We obtain the colu mn density 
measurements from the iSavage et al.l (|2003l) catalog of 
~ 100 sight lin es, which excludes h igh velocity O VI and 
is based on the lWakker et all (|2003T ) data analysis. Their 
four sight lines nearest to our off-cloud sight line are those 
toward PG 1626+554 (log iV = 14.25 ±0.09 ±0.08, where 
the first uncertainty is the combined statistical and con- 
tinuum placement 1 a error in the log of the column 
density; the second is the 1 a systematic error, which, in 
this case, is not necessarily correlated from one observa- 
tion to the next), Mrk 487 (log AT = 14.25 ± 0.11 ± 0.04), 
Mrk 290 (logiV = 14.21 ± 0.12 ± 0.05), and Mrk 876 
(log AT = 14.43 ± 0.02 ± 0.05). The average Milky 
Way O VI column density on these 4 lines of sight is 

i.97tS#18i2 * io 14 cm- 2 . 

The average angular separation between these sight 
lines and our own is 7.3° and the maximum is < 10°. The 
average deviation in the mean column densi ty for pairs of 
sight l ines having ~ 7° of separation in the ISavage et al.l 
(|2003H catalog (their figure 11) is 24%. Following the 
logic of IShelton etaLl (j2007ft . the rms deviation between 
the mean of the parent population and the mean of a 
sample of four observations drawn randomly from the 
parent population should be ypk : /4 times as large as the 
average deviation between 2 sight lines, thus 11%, and 
the intrinsic fluctuation between column densities along 
different directions should be a factor \/tt/2 times the 
average deviation between 2 sight lines, thus 21%. Com- 
bining these terms in quadrature yields 24%. This is 
larger than the average statistical plus continuum place- 
ment error and the average systematic error obtained 
from the four individual column density measurements. 
Here, we add the 3 sources of uncertainty in quadrature 
in order to obtain an estimate of the uncertainty in the 



estimated column density for our direction. The column 
density and combined error are 1.97j^5i x 10 14 cm -2 . 

The column density of local material must be sub- 
tracted from the line of sight column density. The 
ISavage fc Lehnerl tf2006h survey of O VI toward nearby 
white dwarfs includes 9 sight lines that are within 20° of 
our off-cloud direction. Ions of O VI we re not detected on 
4 of t hese sight lines. In those cases, ISavage fc Lehnerl 
(2006J) listed 2 a upper limits. We treat these cases as 
if the la upper limits were half of the 2a upper limits 
and the detected values were 0, so that we can find the 
average for the 9 sight lines. The resulting average vol- 
ume density of O VI is l-St^; 2 x 10~ 8 cm " 3 . Although 
there are reports fi.e,. IBarstow et al"]|2009|) claiming that 
only a few of the previous O VI detections for nearby 
sight lines can be shown to be interstellar, as opposed to 
photospheric or ambiguous (generally due to the similar- 
ities between the observed velocities and the velocities of 
photospheric lines in the spectr a), a culled datase t is not 
yet available. Thus, we use the [Savage fc Lehnerl (2006) 
dataset. 

We estimate the maximum extent of the Local Bubble 
from a survey of the Loc al Cavity's wall, in this case the 
Na I and Ca II survey of I Welsh et al.l (|2010f ). In Figures 
14 and 17 of that paper, the Local Cavity extends 75 pc 
in the direction of our observation. Taking this and the 
average volume density yields an estimate of the Local 
Bubble's column density of O VI of 3.41 2 ; 8 , x 10 12 cm~ 2 . 
This is so small compared with the full column density 
through the Galactic disk and halo that approximations 
are inconsequential. Subtracting this value from the full 
column density yields a halo O VI column density of 
Novi = 1-9 ± 0.5 x 10 14 cm -2 . Again, the unrounded 
measurement values are used in the calculations, but the 
rounded results are reported. 

2.4. The Halo O VII and O VIII Intensities 

2.4.1. Intensity Measurements from XMM-Newton 
Observations 

There are 33 archival XMM-Newton observations 
within 5° of our FUSE pointings that have at least some 
EPIC-MOS exposure. All but 3 of the observations were 
badly affected by soft proton contamination. We have 
processed the data fr om these 3 observ a tions using the 
method described in iHenlev fc Sheltonl (|2010T l . We re- 
moved times affected by soft proton flares, and also times 
when the solar wind proton flux exceeded 2 x 10 s enr 2 
s -1 (the latter step was to reduce contamination from 
solar wind charge exchange X-rays; see Section 2.4.2, be- 
low). We extracted spectra from the blank sky regions 
of the EPIC-MOS chips, and measured the intensities in 
the O VII Ka triplet (569-574 eV) and O VIII Ka line, 
accounting for the effects of residual soft proton contam- 
ination and the extragalactic background (the spectrum 
of which we assumed to be 10.5(.E/keV) _146 ph cm- 2 
s- 1 sr- 1 keV -1 ; iChen et al.lll997D . For the portions of 
observation 0200750601 {£ = 90.0°, b = 38.4°) not af- 
fected by soft proton flares, the solar wind proton flux 
exceeded our standard filtering threshold of 2 x 10 8 cm- 2 
s . However, the proton flux was not unusually large, 
and did not exceed 2.8 x 10 s cm- 2 s _1 . We therefore did 
not carry out proton flux filtering on this observation. 
The X-ray intensity measurements are listed in Table [4] 
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Table 4 

O VII and O VIII Intensities from XMM-Newton Data 



Obs. ID 


Obs. date 


£ 


6 


lovn 
(LU) 


lovin 
(LU) 


0049540401 
0200750601 
0302310101 


2003- 08-19 
2005-09-07 

2004- 09-29 


86.567 
90.020 
86.462 


40.932 
38.406 
41.115 


11 6 +0 - 8 

6-818:2 


y + 0.5 

°- ' -0.4 

1 0+ ' 4 

n 9 + 0.7 
z,z -0.5 



Table 5 

Intrinsic Halo O VII and O VIII Intensities for 
I = 90.0°, 6 = 38.4° direction 



Note. — For these measurements, the extragalactic back- 
ground component was removed. During the estimation of its 
strength, the background was assumed to have been partially ab- 
sorbed by Galactic material, whose H I column density was mea- 
sured by the LAB survey. 

2.4.2. SWCX Abatement 

The measured O VII and O VIII intensities differ from 
one sight line to the next by factors of ~ 2 and ~ 4, 
respectively, in spite of the fact that the ROSAT 3/4 
keV map of this region shows a relatively uniform sur- 
face brightness. The variation found between the three 
XMM-Newton observations is not due to absorbing ma- 
terial along the lines of sight, as the obscuration toward 
I = 90.0°, b = 38.4° is insufficient to explain the rela- 
tive dimness of its O VII measurement. The variation 
is probably due to flares of contamination by solar wind 
charge exchange (SWCX) X-rays. Such X-rays were gen- 
erated after solar wind ions charge exchanged with neu- 
tral gas in the heliosphere and Earth's extended atmo- 
spher e (|Robertson fc Cravens! 120031 : iKoutroumpa et all 
2006). Observational work has shown that SWCX 
cont a mination can vary significantly as a function of 
time (ISnowden. Collier fc Kuntzl I2004t iFuiimoto et al.l 
[20071 IHenlev fc SheltorJ I2008t ICarter fc Sembavl 120101) . 
even d uring periods having m odest solar wind proton 
fluxes (|Henlev fc Sheltonll2010D . The intensity variation 
that we found among the three Draco pointings, which 
were observed during different years (see Table 2]), is 
within the range of the variation seen among multiple 
XMM-Newton observatio ns of the same direction in the 
IHenlev fc SheltorJ (|2010D catalog and is probably due to 
variability in the SWCX X-ray intensity. With this in 
mind, we take the observation having the minimum oxy- 
gen intensities as that which is least affected by time- 
variable SWCX X-rays, and adopt these measurements. 
Nonetheless, this observation, (£ = 90.0°, b = 38.4°, hav- 
ing Iovn ~ 7 LU and Iqviii ~ 1-0 LU) is subject to 
low-level, slowly-varying SWCX. Note that the pointing 
direction is only 1.2° from our FUSE and FIMS/SPEAR 
direction and that the O VII and O VIII intensities are 
not unusual co mpared with those s een along b ~ 30° di- 
rections in the IHenlev fc SheltorJ (|2010D XMM-Newton 
catalog. 

2.4.3. The Intrinsic Halo O VII and O VIII Intensities 

Low-level SWCX contamination is not the only con- 
tribution to the observed spectra. SWCX X-rays, to- 
gether with Local Bubble X-rays, if they exist, form a 
foreground component, If g , which is not subject to in- 
terstellar extinction. The emission from the Galactic 
halo, It, is subject to extinction having an optical depth 
of r and combines with the foreground component as 
lobs = Ifg+Ih exp(— t) to produce the observed intensity, 



Method of Measuring 
Obscuration 


N H 
10 20 cm- 2 


Ih,OVII 

(LU) 


Ih,OVIII 

(LU) 


LAB Survey 

/i00 M m (0.920 MJy/sr) 

E(B-V) (0.0241 mag) 


1.93 
1.48 
1.65 


6.3+i'l 


1.1 ± 1.2 
1.0 ± 1.2 
1.0 ± 1.2 



Note. — The assumed optical depth affects the re- 
moval of the extragalactic background. For consistency, 
the data were reprocessed whenever the assumed optical 
depth was changed. The absorption cros s-sections are from 
Balucinska-Church & McCammon (1992) with an updated 
He cross-section from lYan et alj (119981 ) . calculated using 

Anders & Grcvcssc (1989) abund ances. 

Line 2: the conversion re lation in Snowdcn ct al. ( 20Q3) was 
used to convert from the Schlcgc l et al.l Q199SI 1 Jinn.irr. inten- 
sity to Nh- 

Line 3: the conversion relation from IGiiver fc Ozell J2009h 
was used to convert the Schlcgcl ct al. 
V) to N H . 



(1998) value of E(B- 



Iobs- The foreground component can be estimated from 
shadowing o b serva tions of nearby clouds (see Table 4 in 
Gupta et all (120091) ] and SWCX models (see Table 4 in 
Koutroumpa et all 1)20071 )) for cases not subject to strong 
SWCX flares. Here, we take the foreground component 
to be If g ,ovil = 1-5 ± 1-0 LU and If g ,ovin = ± 1.0 
LU. The optical depth was estimated in three different 
ways: by using the weighted average of the H I col- 
umn densities on ne arby pointings in t he LAB survey 
Kalberla et "aLl[2005h . using relations in ISnowden et all 



,2000) to convert the 1 00 ^m intensity from the DIRBE- 
corrected IRAS maps ()Schlegel et al.lll998D to N H , and 
using relations in IGiiver fc Ozell ()2009D to convert the 
reddening (E(B-V)) obtained from the DIRBE-corrected 
IRAS m aps to Nh- In each case, the absorpt i on cro ss sec- 
tions of iBalucihska-Church fc McCammon ([1992) were 
then used. Table [5] lists the data from these surveys 
and the intrinsic O VII and O VIII intensities calculated 
by solving for Ih in I i> s = If g + 7/ l cxp(— r). Here, we 

LU and 



+ 1.4 
-1.5 



adopt the average values: Ih ovu = 6.4 
Ih,ovin = 1-0 ± 1.2 LU. 

3. DISCUSSION 

The ratio of Iovi to Novi allows us to estimate the 
electron density (n e ) and thermal pressure (pth) of the 
O Vl-rich gas in the halo. The calculat i on me thod 
is the same as th at used in IShelton et al.1 (|2001D and 
IDixon et al.1 ([2001), but is used here with the intrinsic 
intensity and column density of the O VI ions located 
above the disk. It yields n e = 0.011 ± 0.004 K cm" 3 
Pth/k = 65OOI2600 K cm -3 , where k is Boltzmann's con- 
stant. This pressure estimate and that found for the 
southern off-filamen t direction (7000 - 10,000 K cm" 3 ; 
IShelton et al J 1200 7) are the only pressure estimates for 
extraplanar O Vl-rich gas in which the local contribu- 
tions to the observed quantities were subtracted before 
the pressure was calculated. In addition, several esti- 
mates have been made without the subtraction of the rel- 
atively small local co ntributions to the O VI intensity and 
column density (i.e.. IShelton et al.ll200lt P th /k ~ 5300 
to 14,000 K cm" 3 , with the range reflecting the un- 
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certa inty in the location of the H I gas; iDixon et al.l 
[200lt Pfh/k ~ 20, 000 K cm" 3 for T = 10 5 3 K gas; 
ISheltonl[200l P th /k ~ 3700 to 4 900 K cm" 3 , depending 
upon the location of the H I gas; IDixon fc Sankritll2008l : 
Pth/k ~ 5300 to 7400 K cm" 3 for the Fairall 9 direction 
and 3000 to 3600 K cm~ 3 for the NGC 625 direction, 
with the ranges reflecting the uncertainties in the loca- 
tions of the H I gas). 

The thick-disk/halo O Vl-rich gas is expected to reside 
between heights of \z\ ~ 160 pc (the height of the shad- 
owing cloud that was used in the observation of local O VI 
photons and discussed in Subsection l2.2.3l) an d ~ 23 00 pc 
(the scale height of O VI ions. I Savage et al.l (|2003fl ) and 
has a much larger thermal pressure than had been es- 
timated for this region based on observations and the 
assumption of hydrostatic equilibrium {Pth/k < 2200 K 
cm~ 3 ; lFerrierdll998l ). The discrepancy could be due to 
a lack of hydrostatic balance, as would occur if disrup- 
tive events occurred too frequently for the thick disk to 
re-relax. The discrepancy could also be due to the as- 
sumption that the gas is at its CIE temperature. The 
calculated pressure is proportional to temperature di- 
vided by a function of temperature and has a minimum 
of p^/k = 3300±1300K cnr 3 at T = 9.6 x 10 4 K. Alter- 
natively, the discrepancy for our direction can be reduced 
by relaxing our assumption that the obscuring material 
along the sight line resides below the extraplanar O VI. If, 
instead, it lies above, then the extraplanar O Vl's inten- 
sity would be 3480l^oo LU ( from Subsection [2X3]), im- 
plying a thermal pressure oipth/k — 4700jZi|oo K cm ~ 3 - 

The ratios of the intrinsic extraplanar intensities of 
C IV, O VI, O VII, and O VIII are useful tests for halo 
models. Here we use them to test several models, be- 
ginning with a model in which gas is heated to a high 
temperature and then allowed to cool radiatively. This 
model is plausible, given that its predictions have been 
found consistent with the ratios of O VI, O VII, and O VIII 
col umn densities observed for the PKS 2155-304 sight line 
bv iHeckman et all (|2002h . An application of the radia- 
tive cooling model involves g as that was a c creted from 
intergalactic space. Although iHenlev et ail (|2010f l show 
that an extended hot halo of accreted gas (r ~ 10s of kpc) 
may be too faint to account for the X-ray emission seen in 
the XMM-Newton band (including the O VII and O VIII 
lines) , the model X-ray luminosities used bv lHenlev et al.l 
(2010) excluded emission from near the disk. If accreted 
material near the disk contributes to the observed O VII 
and O VIII emission then the model is not ruled out. 
Here we predict theoretical intensity ratios and compare 
w ith our observations. 

iLei. Shelton fc Henlevl ()2009f ) presented a simple 
model in which radiative cooling of hot, shocked inter- 
galactic gas results in a temperature-stratified layer of 
accreted gas. In such a model, the intensities in the C IV, 
O VI, O VII, and O VIII lines can be calculated from 

-(l)(H(if)/^- » 

where rih and n p i are the volume densities of hydrogens 
and positive ions, respectively, s is fo r isochorically 
cooling gas (|Lei. Shelton fe Henl ey 2009) and 1 for iso- 

baricically cooling gas (E dgar fc Chev alier 1986), \ ^ is 



the unknown rate at which particles are accreted per unit 
cross sectional area, e(T) is the emission coefficient, and 
Ajv(T) is the cooling function fo r radia tively cooling gas. 
Following iSutherland fc Dopital (|1993f ). A^ has units of 
ergs cm 3 s _1 and is defined such that the energy loss 
per unit volume per unit time equals n e n p iAj\r(T), where 
n e is the volum e density of electrons . We take An(T) 
from Table 6 in ISutherland fc Dopital (|1993f ). which ap- 
proximates the ionization levels of the gas as those in 
CIE p lasma. Solar abundances from I Anders fc Grevessd 
(1989) are used. We take the emission coefficients from 
the Xspec databases for the Raymond & Smith and 
APEC models for radiation from optically thin thermal 
plasmas in collisional ionizational equilibrium. The ele- 
mental abundances , which are assumed to be solar from 
lAnders fc Grevessel (|1989f ). and the fraction of atoms in 
the relevant ionization stage are factored into these coef- 
ficients. The emission coefficients have units of ergs cm 3 
s~ 4 and are defined such that the energy loss per unit 
volume per unit time is e(T)n e n#. Because An{T) and 
e(T) are tabulated at discrete values of T, we approxi- 
mate the intensity integral with a sum. 

Not knowing the value of -^ ^ , we report only the ra- 
tios of the intensities for comparisons with the observa- 
tions. The model ratios of Iqiv '■ lovi '■ lovn '■ Ioviii 
are 0.13 : 1.00 : 0.055 :> 0.26 from the Raymond & 
Smith coefficients and > 0.068 : 1.00 : 0.057 :> 0.26 
from the APEC coefficients. Upper limit signs appear 
for ions whose emissivity functions extend to higher or 
lower temperatures than were covered in the tables. (The 
Sutherland & Dopita tables run from T = 10 4 to 10 8 5 K, 
while the Raymond & Smith tables extend from T = 10 4 
to 10 8 K and the APEC tables extend from T = 10 5 to 
10 8 ' 9 K.) These ratios were calculated from intensities in 
units of photons sec -1 cm -2 sr , rather than intensi- 
ties in units of ergs sec -1 cm~ 2 sr -1 . The derived ratios 
have far too little C IV and O VI flux, or conversely, 
far too much O VII and O VIII flux for consistency with 
the observationally derived ratios 1.4 ± 0.6 : 1.0 : 1.3 ± 
0.5 x 10~ 3 : 2.2tf e 5 x 10~ 4 a problem that also plagues 
the extraplanar gas near the filament in the southern 
galactic hemisphere (Iciv ■ lovi ■ lovn ■ Ioviii = 
1.0 ±0.4 : 1.00 : 1.3± 0.2 x 10~ 3 : 3.4±g;g x 1Q- 4 , fr om in- 
tensities presented in ILei. Shelton fc Henley 2009). This 
problem cannot be remedied by adjusting the gas phase 
metal abundances, although it can be reduced by de- 
creasing the temperature of the upper limit of integra- 
tion in the intensity integral. The temperature of the 
freshly accreted gas should be around 10 6 5 K before 
it begins to cool (|Lei. Shelton fc Henlevl 120091 ) . Setting 
the upper limit of integration to T = 10 6 ' 5 K results in 
Iciv ■ lovi '■ lovn ■ Ioviii = 0.13 : 1.00 : 0.037 : 0.031 
and > 0.068 : 1.00 : 0.041 : 0.027 for the Raymond & 
Smith and the APEC coefficients, respectively. The pre- 
dicted ratios of lovn to lovi and Ioviii to lovi are 
still far too large. Our model could be tuned by fur- 
ther adjusting the upper temperature cut off; however, it 
would be more logical to interpret this model as a source 
of some of the background O VII and O VIII, while the 
majority of the C IV and O VI on the off-cloud sight line 
derive from other phenomena. 

Localized, sporadic injections of energy preferentially 
produce C IV and O VI and could explain the uneven 



8 



Shelton et al. 



distributio n of O VI ions acro s s the sky (see O V I distri- 
butions in ISavage et all 120031: iDixon et all I2006D . Sim- 
ulations of impulsive injections of energy in the form of 
supernova explosions have been done for explosion en- 
ergies of 5 x 10 50 ergs occurring in an ambient medium 
that has a combi ned magnetic and cosmic ray pressure 
of 1800 K cm" 3 (jSheltonl 1200(1 . These large injections 
of energy engender shocks that heat and ionize the gas. 
Although the early post-shock temperatures are above 
3 x 10 6 K, the expanding spherical shockfronts weaken 
with time, reducing the post-shock temperature. As this 
happens, the gas inside the bubble cools by adiabatic 
expansion. Later, it cools primarily by radiative cool- 
ing. Each of these processes skews the temperature dis- 
tribution in favor of C IV and O Vl-rich gas rather than 

VII and O VHI-rich gas. When averaged over time 
and space, the resulting bubbles have ratios of median 
intensities of Iciv ■ lovi '■ lovn ■ Iqvttt = 0-39 : 1.0 : 
3.4 x 10~ 4 : 3.7 x 10~ 5 (assuming lAnders fe Grevessg 

1 !)>■>!) abundances). These ratios were calculated from 
a Monte Carlo simulation that included the luminosity 
of individual SNR bubbles as a function of time, varia- 
tion in ISM conditions and therefore SNR characteristics 
as a function of height above the plane, and the likeli- 
hood that a sight line would pass though any given SNR, 
based the SNR's size, height above the midplane, and 
the height distribution of progenitor stars. They, there- 
fore, represent the all high-latitude-sky average, which is 
serendipitous as our observed intensities are also near the 
average of those observed across the high-latitude sky. In 
contrast with the accreting gas model, the impulsive in- 
jection model predicts too little emission from O VII and 
O VIII per photon emitted by O VI. 

Another possible source of extraplanar O VI, O VII, 
and O VIII is hot gas that has welled up from the 
Galactic disk. In the h ydrodynamical simulations of 
Uoung fc Mac Low! (|2006f ). SN explosions stir the inter- 
stellar medium and drive hot gas into the thick disk. 
These simulations, which modeled 80 Myr of ISM evolu - 
tion, were extended to 155 Myr for lHenlev et al.1 (|2010H . 
(Note that these times include the time taken for the 
initial conditions to be eradicated.) From them and the 
Raymond & Smith spectral code, iHenlev et al.l ()2010f ) 
created synthetic spectra. The synthetic spectra from 
time periods after 140 Myr were found to be at least 
as bright in 0.4 to 2.0 keV photons as the spectra ex- 
tracted from XMM-Newton observations after removing 
the local contribution. Thus, such upwelled hot gas is 
a probable source of O VII and O VIII. However, model 
Iciv ■ lovi '■ lovn ■ Ioviii ratios are not yet available 
for comparison with our observations. 

Furthermore, there have been reports of en- 
hanced X-ray inte nsities associated with high veloc- 
ity clouds, HVCs mirth et al.l fl985t IKerp et all llMi 
iBregman et al.l 12009ft . Our line of sight passes near the 
intermediate and high velocity material associated with 
the Draco Nebula and Complex C, crosses some fast mov- 
ing material (see Section l2.2.2[) . and hypothetically could 
encounter materal left behind by an HVC, given that sim - 
ulated HVCs shed material (jHeitsch fc Putmanl [2009h . 
Although intensity predictions from simulations are not 
yet available, we can assume that several known physical 
processes may be active. Supersonic collisions between 
high velocity clouds and thick disk/halo gas would shock 



heat the gas. The clouds may experience turbulent mix- 
ing and evaporation if they are passing through a hot 
ambient medium. Radiative cooling will also operate. 
iSlavin et al.1 (|1993| ) analytically calculated the C IV and 
O VI intensities resulting from model turbulent mixing 
layers. The Iciv ■ lovi ratios from their predictions 
ranged from 1.7 to 16, depending upon the assumed con- 
ditions. The low end of the predicted Iciv ■ lovi ratio 
is within the la error bars of the observationally deter- 
mined ratio and is from models having m i xed g as tem- 
peratures of T = 10 5 K. iKwak fc Sheltonl (|2010( i hydro- 
dynamically simulated mixing layers and calculated the 
quantities of C IV, N V, and O VI, ions in a time de- 
pendent manner (i.e., not assuming CIE). Subsequently, 
they calculated the ratio of the C IV photon intensity to 
the O VI photon intensity for their primary model, Model 
A, and for a similar model with 10 times as much thermal 
pressure after 50 Myr of evolution, finding it to be 2.4 
and 4.9, respectively (K. Kwak, personal communcation, 
2010). Not only are the ratios of C IV to O VI intensity 
larger in turbulent mixing layer models than in radiative 
cooling and supernova remnant models, but the ratios of 
C IV to O VI column densities are also larger (ISlavin et al.l 
fl^lBs^fivel et alll2006l: IKwak fc Sheltonll2010D . None 
of these papers predicts the intensities of O VII and O VIII 
and such values would depend strongly upon the assumed 
depth of the hot reservoir. However, we can assume that 
the greater the number of mixing layers, the greater the 
lovi '■ lovn and lovi '■ Ioviii ratios will be. 

In the foregoing analyses, we have ignored scattering 
of UV background photons by interstellar C IV and O VI 
ions be cause we expect the scattered intensities to be 
small. iShelton et al.l ()2001f ) estimated the intensity of 
photons scattered by O VI ions, finding it to be much less 
than (i.e., 6% of) the intensity produced by line emission 
for their Case 3b, which is similar to our case in that 
it assumed that most of the O VI photons originated in 
the halo and that there is some line-of-sight obscuration. 
By similar logic, we expect that scattering by O VI ions 
has negligibly affected our observations. To examine the 
magnitude of scattering by C IV ions in co mparison with 
that by O VI ions, we use equation 14 from IShelton et all 



(|200lD for the scattered intensity, I S: 



-i2b j 
A scl 



F(t 



The background continuum specific intensity, I SC i, near 
1550 A is 4.1 x 10~ 21 ergs cm -2 s _1 sr _1 Hz -1 ()Bowverl 
119911 ). which is only 1/10 of the intensity near 1030 A, 
4.4 x 10" 20 ergs cm " 2 s" 1 sr" 1 H z " 1 dMathis et al.lll983l 
and assumed in the IShelton et al.l (|2001[ ) analysis! . If the 
line profile is thermally broadened and the gas is in CIE, 
then the velocity spread parameter, b, of C IV is 2/3 of 
that of O VI. Its wavelength, A, is obviously larger. We 
expect the column density of C IV ions to be around 
60% of tha t of O VI ions, giv e n the ensemble statistics 
reported in ISavage fc W akker ( 2009). Thus, the optical 
depth of scatterers, r , will be smaller, and F, which is 
a function of r Q , will be smaller for C IV than for O VI. 
The variable a is related to absorption by intervening 
dust and is 6% larger at 1550 A than at 1030 A. As a 
result of these factors, the intensity scattered by C IV 
ions is even smaller than that scattered by O VI ions. 
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In this paper, we discuss measurements of the C IV, 
O VI, O VII, and O VIII intensities. The O VI, O VII, 
and O VIII data are ne wly presented here, while the C IV 
intensity is taken from iPark et all (j2009| ) . Each of the 
measured intensities was taken from I = 90°, b ~ 40°, a 
region that is near to but not toward the Draco cloud, 
and is fairly typical of high latitude values. 

We have subtracted the contributions made by local 
diffuse gas and compensated for the absorption loss due 
to intervening material, in order to obtain the intrinsic 
intensities of photons that originated above the plane. 
While this region is often called the halo, especially when 
discussing the diffuse X-ray background, most of the ob- 
served ions reside in the thick disk. We find a high pres- 
sure pth/k = 6500±26o° K cm ~ 3 > see Section [3] for the 
O Vl-rich gas in this region. It is consistent with the 
measurement of extraplanar O VI made in the southern 
Galactic hemisphere and is greater than expected from 
static models of the thick disk. Such a high pressure 
could result from the thick disk being disturbed on a 
faster timescale than the relaxation timescale. In ad- 
dition, if the true temperature of the gas is much less 
than the CIE temperature, then the true pressure may 
be somewhat less than that derived from our measure- 
ments. 

Each of our ions traces hot gas, although they sample 
different, sometimes overlapping temperature subregimes 
(C IV traces ~ 1 x 10 5 K gas, O VI traces ~ 3 x 10 5 K gas, 
O VII traces 1 x 10 6 K gas, and O VIII traces — 3 x 10 6 K 
gas if the gas is in collisional ionizational equilibrium). 
The hot gas in these subregimes is thought to be causally 
and sometimes physically related. One relationship oc- 
curs through radiative cooling, because gas that was 
heated to a few million degrees will cool through the 
other subregimes if left undisturbed. Another relation- 
ship occurs through energy transfer, either by thermal 
conduction or by mixing. If very hot gas abuts signif- 
icantly cooler gas, then a C IV and/or O Vl-rich zone 
develops between the hotter (O VII and/or O VHI-rich) 
gas and the cooler gas. 

The ratio of intensities provides a diagnostic of the 
heating mechanism. Mechanisms that heat gas to several 
times 10 5 K favor C IV and O VI, while mechanisms that 
heat the gas to higher temperatures boost the relative 
O VII and O VIII intensities. It should be noted that 
atomic physics favors the resonance lines of lithium-like 
ions (i.e., O VI 1032, 1038 A and C IV 1548, 1551 A 
emission) over the K a lines of helium-like and hydrogen- 
like ions (i.e., the O VII "triplet" between 561 and 574 eV 
and the O VIII 653 eV line), if all other factors, such as 
abundance, are equivalent. Therefore, when we indicate 
that heating the gas to very high temperatures favors 
O VII and O VIII emission, we mean that Iovn /lovi an d 
Ioviii / lovi are greater than in other models. We do 
not mean to imply that our O VII and O VIII intensities 
are brighter than our O VI intensity, because they are 
not. 

The O VI intensity is the reference point for all of 
our comparisons. Our intrinsic, halo O VI intensity 
is 4770±^[J LU (see Subsection |2"X51) and the I C iv ■ 
lovi '■ Iovn '■ Ioviii ratios are 1.4 ± 0.6 : 1.0 : 
1.3 ± 0.5 x 10~ 3 : 2.2±H x 10~ 4 (see Section©. For 
the moment, we ignore C IV, which we will return to 



later. The lovi '■ Iovn ■ Ioviii ratios lie between the 
extremes formed by two models, one of which is rela- 
tively preferential to O VII and O VIII while the other 
is relatively preferential to C IV and O VI. The first 
model is fairly simple. It heats gas to a high temper- 
ature (10 6 5 K), then allows it to cool radiatively. This 
model, which may correspond to halo gas resulting from 
a smooth accretion of intergalactic gas, produces ratios 
of 0.13 : 1.00 : 0.037 : 0.031 (assuming Raymond & 
Smith coefficients; the ratios found from APEC coeffi- 
cients are also presented in Section©. The second model 
was made from simulations of the bubbles blown by spo- 
radic and isolated injections of energy. The particular 
model employed was for bubbles blown by supernova ex- 
plosions. The model ratios take into account the differ- 
ent numbers of photons produced at different stages of 
the bubble's life, including the short-lived early period 
when the shock-front heats the gas to temperatures well 
above 3 x 10 6 K, the longer-lived stage when the shock- 
front heats the gas to a few times 10 5 K, and the final 
much longer period when the SNR bubble radiatively 
cools. Thus, this model includes the effects of inject- 
ing energy to create a wide range of gas temperatures, 
in addition to radiative cooling. It produces ratios of 
0.39 : 1.0 : 3.4 x 10~ 4 : 3.7 x 10" 5 (see Section©. The 
observed lovi '■ Iovn '■ Ioviii ratios are between the 
two comparison models. 

We can conclude that in addition to some very hot gas 
(Icze~3x 10 6 K), a broad temperature profile is needed. 
Not only could a broad temperature profile result from 
supernova bubbles, but it may also result from mixing 
or thermal conduction between very hot gas and cooler 
gas, injections of moderately hot fountain gas into the 
"halo" , or other phenomena. 

Interestingly, neither model produced more C IV pho- 
tons than O VI photons, yet the observations show that 
C IV is at least as intense as O VI in both the Draco re- 
gion and the southern filament region (see iShelton et al.l 
l2007f )Fl Adjusting the C:0 abundance ratio cannot solve 
this problem, as it woul d require raising it from 0.4 
(| Anders &; Grevessel 1 19891 ) to 1.6, which is unrealistic. 
Instead, the large Iciv '■ lovi ratio indicates that there 
exists a strong mechanism for producing C IV photons. 
Turbulent mixing between hot and cool gas is one such 
mechanism. 
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